Correlations in atomic systems: Diagnosing coherent superpositions 



O 

o 

(N 



(N 



o 



o 
o 



> 
m 
vn 
in 
m 
o ^ 

^ : 
o . 

a ■ 



I 

O 

o 



X 
J3 



Radka Bach and Kazimierz RzazewskQ 
Center for Theoretical Physics, Polish Academy of Sciences, al. Lotnikow 32/46, 02-668 Warsaw, Poland 

(Dated: February 2, 2008) 

While investigating quantum correlations in atomic systems, we note that single measurements 
contain information about these correlations. Using a simple model of measurement — analogous to 
the one used in quantum optics — we show how to extract higher order correlation functions from 
individual "photographs" of the atomic sample. As a possible application we apply the method to 
detect a subtle phase coherence in mesoscopic superpositions. 



Correlation functions provide a useful tool for investi- 
gating complicated states, both classical and quantum. 
For example, measuring intensity-intensity correlations 
for photons has been a significant development in astro- 
physics and investigation of the second order correla- 
tion function has led to understanding of antibunching of 
photons 2]. For atoms, correlation functions might prove 
to be as useful as they have been for photons. There are 
certain important differences, however. For light, if the 
underlying quantum state is unknown, the knowledge of 
a correlation function of a given order does not determine 
lower order correlation functions and a separate experi- 
ment must be set up to measure the latter j3, l^] . This 
is not true in atomic case. For atoms a superselection 
rule holds, which prohibits the atomic state to exist in 
a superposition of Fock states with different number of 
atoms. This proves to be a very strong condition and 
implies, as will be shown in more detail, that there exists 
a unique hierarchy of correlation functions. An impor- 
tant consequence is that it is possible to extract informa- 
tion about higher order correlation functions, not only 
density, from typical measurements performed in atom 
optics. In this Letter we derive this hierarchy of cor- 
relation functions and provide a simple prescription for 
extracting correlation functions of interest from typical 
measurements, which we model in an analogous way to 
quantum optics. As an illustration we use the method as 
a diagnostic tool to distuinguish a coherent superposition 
against a simple mixture. 

There exists a number of possible detection schemes in 
atom optics - via resonance fluorescence, non-resonant 
imaging or ionizing the atoms and detecting the result- 
ing electron current, just to name a few - and each of 
them requires a slightly different theoretical approach . 
In particular one cannot claim that all measurements of 
atomic samples correspond to norniarly ordered corre- 
lation functions and in that respect atom optics differs 
from its photonic counterpart, in which all correlation 
functions are normally ordered since light is detected via 
absorption. On the other hand, however, the notion of 
coherence, that has been originally introduced for pho- 
tons in the context of normally ordered correlation func- 
tions and their factorization jy| , is a very useful tool and 
it would be desirable to classify atomic systems in an 



analogous manner. Therefore we will further assume that 
the measurement of the atomic cloud is done in such a 
way that normally ordered correlation functions are rele- 
vant, even though it might be an idealization of an actual 
detection scheme. 

Consider a state composed of Af atoms and described 
by a ket {tp) for a system with a second-quantized field 
\t'(x). The diagonal part of the r-th order correlation 
function is then defined as: 

i, |*t(xi) . . . ¥{^r)'^{^r) ■ ■ ■ *(xi)| V^) (1) 

and, within the assumptions described above, is propor- 
tional to the joint probability of detecting r atoms: one 
at position xi, one at X2 and so on (of course the po- 
sitions need not be all different). Time dependence has 
been dropped for simplicity: we assume that all atoms 
are measured at the same instant. In other words we 
neglect the state's dynamics during the time of measure- 
ment. 

Provided the system is in a Fock state, i.e. is 
an eigenstate of the total number of atoms operator 
J\f = J dx 5'^(x)^(x), a simple algebra shows that the 
following relation holds both for bosonic and fermionic 
systems: 
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Thus, a correlation function of a given order contains all 
lower order correlation functions. Such a hierarchy exists 
for any given state, but only for Fock state is the coeffi- 
cient of proportionality so simple and general. Assump- 
tion that the state is in a Fock state is valid for most 
experiments performed with cold atoms, Bose-Einstein 
condensates in particular, since they are almost perfectly 
isolated from the environment. Moreover, if one was to 
allow for interactions with the particle reservoir, the re- 
sulting state is a mixture of Fock states with different 
number of atoms: p = X^J^o '""1^") (^"1 ^^'^ exten- 
sion of Eq. [21is straightforward: 

dxr G^ixi, . . . ,Xr) = 
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(n-r + l)G[;-i)(xi,. 



-i) (3) 



density cfiUx) = Ql',l^'{x), since the first one also car- 
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Other states, allowing for coherence between states with 
different number of atoms, violate the superselection law 
that stems from the baryonic charge conservation j2i] and 
thus are unphysical for atomic systems. 

Joint probabilities of detecting atoms at xi , X2 , . . . , 
form hierarchy analogous to Eq. Inland therefore: 

p[;;(xi, . . . ,x.) = i:^(^G;;;(xi, . . . ,x.) (4) 

Let us now investigate the experiment in greater de- 
tail and concentrate on an ideal case first, i.e. when the 
detectors have infinite spatial resolution and are able to 
detect individual atoms. We will also assume that the 
detectors' mesh is so tight that each and every atom gets 
detected. Then, photographing the state produces 
a collection of A/" positions £,i,£,2, ■ ■ ■ at which atoms 
were found at the i-th shot. Repeating the measurement 
W times follows the corresponding probability density 
and therefore yields a discrete estimate of the Af-th cor- 
relation function: 



,(AA,W) 



(xi,...,x^) = -^AA! X{5{^,~i 



(5) 



i=l 



Of course in the limit of repeating the experiment in- 
finitely many times the original correlation function is 
recovered: 



w- 



, y^^f) (6) 



To obtain the density out of these measurements one 
needs only to integrate the constructed correlation func- 
tion over M ~ 1 coordinates according to the Eq. [3 (It 
is assumed that this equation holds true also for the dis- 
crete estimates of correlation functions.) Since it doesn't 
matter over which M —1 variables the integration is per- 
formed one obtains a very intuitive result: 
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(7) 



ries information about correlations in the system [Sj. 

Second order correlation function can be extracted 
from the measurements' outcome in an analogous way 
via integrating Q^^^'^^ (xi , . . . , x^r) over Af—2 variables. 
After exploiting the symmetry trick, just as before, we 
obtain: 
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EE^(^-^)'^(^'-^')(i-'5,,,o (8) 

3 = 1 j' = l 

This idea has been already used in in which the au- 
thors propose to look at spatial noise correlations in or- 
der to extract second order correlation function. Similar 
method can be applied to extract even higher order cor- 
relation functions. 

To describe realistic experiments — with inevitably 
imperfect detectors — we will incorporate the following 
strategy: For each outcome of a realistic measurement, 
the set of all possible ideal realizations (discussed above) 
that are identical from the point of view of imperfect 
detectors, is found; then, for each of these realizations the 
correlation function of interest is calculated and finally 
averaging over these realizations takes place. 

The outcome of the i-th experiment is modelled by a 
set of V numbers denoting how many atoms were found 
by each of the detectors: S^, S2, . . . , 5^. Here, the de- 
tectors are assumed to have a finite spatial resolution 
(each of them covering a volume V) and to fill all avail- 
able physical space. Furthermore, they are character- 
ized by the conditional probability p(r]\S) of the pres- 
ence of 77 atoms while 5 were actually detected. De- 
noting: »7*'(S) := J2riPi^\^)''l'^ (^°^ example, if each 
atom is independently detected with a probability p then 

^ = and ^ = + i^±fi^), first- 

and second-order correlation functions take the form: 

w V 

e^^'^'(x)-^EE^^-W' (9) 
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The density is thus just a sum of 5-spikes at positions 
at which atoms were detected in all experiments. Note 
the self averaging phenomenon: typically the number 
of atoms A/" is of the order of hundreds of thousands 
and one does not need to repeat the experiments very 
many times to obtain a reasonable density approxima- 
tion. However, in general it is not sufficient to perform 
the experiment only once, no matter how many atoms 
there are. The "density" obtained in a single measure- 
ment, Q'^^'^^\x), is -apart from the systems exhibiting 
coherence up to the Af-th order- different from a true 
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A,(x') (rf(E^)-r^(E^)) + J] A,,(x')r;(SJ,; 



Here, Aj (x) = V^^ if x is within the j-th detector and 
zero otherwise. Note that the knowledge of the under- 
lying quantum state is not required here and still lower 
order correlation functions can be extracted. 
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One more remark about experimental reality. A typi- 
cal measurement performed with cold atoms involves so 
called column averaging, in which information about cor- 
relations in the direction along the propagation of the 
illuminating light pulse is lost. The scheme of extract- 
ing correlation functions presented above is not sensitive 
to whether the original high-order correlation function is 
column averaged or not, and in this sense this process is 
not an obstacle. But if one starts from a column averaged 
high-order correlation function one ends up with column 
averaged low-order correlation functions. In many situ- 
ations this is not satisfactory, however, because such an 
averaging kills interesting phenomena. Consider, for ex- 
ample, the recent experiment of collapses and revivals 
of matter wave field of a Bose-Einstein condensate pTij |. 
Only because of the column averaging were the collapses 
of the interference pattern observed in a single shot. If 
one repeated the experiment with 2-dimensional lattices, 
interference patterns would be observed at each mea- 
suerement no matter what the free evolution time was 
and they would vanish only when averaged over many 
realizations. 

Let us illustrate this general method of extracting 
higher order correlation functions with an example of a 
coherent superposition of states composed of J\f atoms of 
the form: 



(a) 



1 



ICOH 



V2 



{\Af,Q) + \QM)) 



(11) 



Such a state provides a very important link between 
quantum and classical worlds, is essential in the theory 
of quantum information, and is currently being a subject 
of intensive research . We will assume that the wave- 
functions associated with the two modes are the ground 
states of spatially separated potential wells: fL{x) and 
ipfi{x). The r-th order correlation function is: (a) if 
r < TV: 

G'coff (■''1' -^21 ■ • ■ J ^r) = 



2{Af~r)\ 
(b) if r = Af: 



1[\^l{x.o\' +i[\vR{x.)n , (12) 



\i=l 
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(xi, . 



,XAr) = — 



M AT 



(13) 

and (c) if r > A/" the correlation functions vahish inden- 
tically. 

Explicit formulas for correlation functions allow us to 
simulate the experiment easily. Since each measurement 
produces a set of numbers drawn according to the corre- 
sponding joint probability function, so can it be done in 
the computer with the help of random number generator. 
This is achieved step by step: the first atom is drawn with 
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FIG. 1: Low-order correlation functions for a superposition 
state \i^)coH with M = 100 atoms: (a) density obtained in a 
single measurement, (b) density averaged over 50 shots and 
(c) second order correlation function averaged over 50 mea- 
surements. The wavefunctions 'Pl,r{x) were initially Gaus- 
sian wavepackets centered at = ±q: [a denotes the unit 
of length) with width a = 0.035q to make them practically 
orthogonal, but they were allowed to expand freely for time 
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t = 0.01^^^ before the detection process. The solid line at 
(a) and (b) shows the expected density profile, G^qqu{x). 



the probability density given by the second - 

with the probability density GcQff{xi — £_i,X2), where 
for xi a number obtained in the previous step is substi- 
tuted and so on 12J . In this way the full set , ^2 , ■ • ■ , £.j\f 
is obtained. 

Figure n shows density and second order correla- 
tion function obtained via the hereinbefore prescription 
from the numerically simulated measurements. Fig.^a) 
shows a typical result of a single measurement, in which 
all atoms are detected in one of the wells. Averaging over 
many experiments leads to the correct density profile, as 
shown at Fig.^b). Fig.QJc) shows the second order cor- 
relation function constructed from single measurements 
and then averaged over many shots. 

It is instructive to consider another problem: suppose 
we want to distinguish between the coherent superposi- 
tion defined above and the mixed state: 



(|AA,0)(AA,0| + |0,AA>(0,AA|) 



(14) 



(Similar problem for four ions in a trap was experimen- 
tally investigated in 13].) These states have all correla- 
tion functions identical apart from the one of the order of 
the number of particles: for mixed state the Eg. 1 121 holds 
even when r = J\f. It follows immediately that in order to 
distinguish between these states an experiment measur- 
ing all atoms must be performed in conditions in which 
the spatial wavefunctions overlap, otherwise the interfer- 
ence term in Eq. 1131 will vanish. Since the correlation 
function of interest lives in the d x A/'-dimensional space 



(d is the physical dimension of the space) , it is practically 
impossible to gather satisfactory statistics of experiments 
to probe this space even for moderate number of atoms. 
Fortunately, there is a way around it. Our purpose is to 
differentiate between two kinds of states and any quan- 
tity that does so is enough. According to what has been 
said, such a quantity must originate from the Af-th order 
correlation function, but it does not need to be the corre- 
lation function itself. Simple calculations show that the 
density of the center-of-mass is indeed such a quantity 
and a direct proof obtained from numerical simulations 
is shown in Figure |21 Note that all atoms must necessar- 
ily be detected to obtain a center-of-mass density, which 
seems impractical at first. However, when performing 
this kind of experiment a very precise control of the in- 
put state is required, in particular the knowledge of the 
total number of atoms, and thus it is natural to discard 
all measurements that detected less than JV atoms and 
this way to overcome finite detector efficiency. 

This method of diagnosing the superposition has of 
course its limitations. When increasing the number of 
atoms, the visibility of fringes in the center-of-mass den- 
sity decreases, as well as the separation between fringes 
and consequently better and better detectors are re- 
quired. Moreover, the sensitivity to the phase fluctu- 
ations from shot to shot increases: if we assume that 
the initial state is proportional to \J\f, 0) -I- exp{i(j))\0,Af), 
where (jj varies from experiment to experiment, it is re- 
quired for the method to work that the phase cj) does not 
change more than over t: / [2M). The number of experi- 
ments needed to distinguish between a coherent superpo- 
sition and a mixed state grows rapidly with the number of 
atoms, but still yields a reasonable number of repetitions: 
under the condition that the peaks' height signifying the 
coherent superposition is 10 times the value of surround- 
ing background, 500 runs are required for M — 2Q, while 
5000 for M = 100 - and bigger coherent superpositions 
are not only difficult to diagnose, but even more, to pre- 
pare. 

Of course the mesoscopic superposition prepared even 
in the best experiments cannot be expected to be ex- 
actly of the form given by Eq. 1111 but rather of the form 
^j,^Q c/cj/c, A/" — fc), where the coefficients Cfe are peaked 
around fc = and k = M . The state we have discussed 
above has been chosen for its simplicity, but the analysis 
may easily be extended to more realistic states. 

Summarizing: the fact that atomic systems cannot ex- 
ist in a superposition of Fock states with different number 
of atoms leads to a very simple hierarchy of correlation 
functions. Basing on that we have provided a prescrip- 
tion how to extract higher order correlation functions 
from typical measurements performed in atom optics. 
We have discussed both the idealized and realistic ver- 
sion of experiments and illustrated the method by ap- 
plying it to a superposition of atomic states. We have 
also proposed a method of diagnosing such a superpo- 




FIG. 2; Absolute value of Fourier transforms of histograms 
of the center-of-mass variable for (a) coherent superposition 
and (b) mixed state for N = 2Q atoms and 1 000 measure- 
ments (histograms are shown in the insets). The spikes at 
k = 7tvt—^^-tIvi , which for our set of parameters corre- 
spond to fc = ±4.24, are fingerprints of the coherence between 
states |7V, 0) and |0,A'>. 



sition (versus an incoherent mixture) via measuring the 
center-of-mass density profile. 
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